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We study measurable quantities of bosonic and fermionic mixtures on a one-dimensional ring. These 
few-body ensembles consist of majority atoms obeying certain statistics (Fermi or Bose) and an 
impurity atom in a different hyperfine state. The repulsive interactions between majority-impurity 
and majority-majority are varied from weak to strong. We show that the majority-impurity repulsion 
is mainly responsible for the loss of coherence in the strongly interacting regime. The momentum 
distribution follows the C/p"* universal behaviour for the high momentum tail, but the contact 
C is strongly dependent on the strength of the majority-impurity and in a different way on the 
majority-majority interactions. The static structure factor of the majority atoms exposes a low- 
momentum peak for strong majority-impurity repulsion, which is attributed to an effective attraction 
not expected for purely repulsive forces. 

PACS numbers: 67.85.-d, 02.30.Ik, 03.75.Hh 


INTRODUCTION 


Non-local measurable quantities of cold atom systems 
play a major role in the experimental probe of their prop¬ 
erties. From the momentum distributions in the first real¬ 
izations of Bose-Einstein condensates 0,i and the time- 
of flight images of phase transitions in optical lattices Q , 
to the probing of the Tonks-Girardeau (TG) gas (an in¬ 
finitely repulsive one-dimensional bosonic ensemble with 
fermionic properties 0,0); experimental techniques have 
given rise to an unprecedented degree of control, manip¬ 
ulation and measurement of atomic systems 0^. For 
one-dimensional systems, theoretical research has devel¬ 
oped analytical and numerical methods to study proper¬ 
ties of those quantities, like the universal asymptotic 
behaviour of the momentum distribution |9l-llll| , which is 
governed by the contact C. This is a notion that cap¬ 
tures all universal properties of such systems even close to 
phase transitio ns |l2ll. and has been recently measured in 
ultracold gases |I3l4l5l |. The structure factor, another im¬ 
portant concept originated in solid-state physics to probe 
crystalline lattices, has recently been studied |I6l4l8l | and 
observed for ID Bose gases via Bragg spectroscopy [l9| . 


More recent advances in experiments deal with mix¬ 
tures of Bose-Bose, Bose-Fermi and Fermi-Fermi cold 
gases 23- Hi ■ It became also possible to realize few-body 
ensembles of bosons 2J, [ 23 | and fermions in different hy¬ 
perfine states (H I s-nd measure the effects of an impurity 
by increasing the number of fermions one by one [27j| . For 
composite mixtures of bosons and fermions several local 
quantities have been studied and have exhibited various 
phases not present in pure ensembles , but some 

advances are still to be done on the side of non-local cor¬ 
relators [31 1. 


In this work we investigate, using a Jastrow-type 
ansatz, some measurable non-local quantities for a few- 


body ensemble of bosons or fermions in a one-dimensional 
ring in the presence of an impurity of the same mass but 
in a different hyperfine state. In these systems, repul¬ 
sive interactions between impurity-majority or - for the 
bosonic case - majority-majority pairs can be tuned via 
Feshbach resonances. For the integrable system consti¬ 
tuted by an impurity in a Fermi sea, the ground state 
wave functions and energies have been exactly obtained 
in [Hj and results for correlations are discussed in [H - 
For a system of bosons, if all interaction strengths are 
equal, the Lieb-Liniger integrable model can be applied 
j34| . However, when the interaction strengths differ and 
the system becomes non-integrable [sH - fH ; '''^e show that 
all measurable non-local correlation quantities exhibit a 
behaviour substantially different from that of the inte¬ 
grable one. In particular, the reduced one-body density 
matrix of the majority atoms, a correlation function that 
shows the degree of coherence in its off-diagonal terms, 
is strongly dependent on the impurity-majority coupling 
strength. At the same time, the high momentum tail 
of the momentum distribution obeys a universal 

behaviour for all cases studied here. The numerical val¬ 
ues for the contact C depend not only on the impurity- 
majority and majority-majority couplings, but also on 
the nature of the atom being considered. For different 
interactions, the static structure factor deviates from the 
integrable case, above and below the phononic behaviour 
of the Tonks-Girardeau gas; most surprisingly, a pro¬ 
nounced peak arises for strong majority-impurity repul¬ 
sion, which is an indicator of a quasi-crystalline structure 
and of effective attractive correlations in a purely repul¬ 
sively interacting ensemble. The results obtained here 
for the correlations via a simple Jastrow ansatz are com¬ 
pared with the exact solutions in the corresponding cases 
where those solutions exist, showing in general very good 
agreement. 
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FIG. 1. (Color online) Schematic depiction of the physical 
system for representative cases: (a) strongly repulsive inte- 
grable case, in which the separation between the three atoms 
is at maximum, (b) The interaction between the majority 
pair is strong, but the impurity-majority interaction is weak; 
therefore, the impurity atom appears “delocalised” across the 
ring, (c) The impurity strongly repels the majority pair but 
the majority-majority interaction is weak, so these atoms tend 
to bunch together on the other side of the ring. 

SYSTEM HAMILTONIAN AND ANSATZ 

We consider a system of 3 atoms on a one-dimensional 
ring of length L with contact interactions and periodic 
boundary conditions. The Hamiltonian reads 



where the lengths are in units of L and energies in units 
of fi^/mL^. The interaction strength gij may be differ¬ 
ent for each pair of atoms i,j and is controllable via 
Feshbach or confinement induced resonances since g = 

with 9.D = 

where m is the atom mass, ^ is the Riemann zeta func¬ 
tion, asD is the 3D s-wave scattering length and a_L the 
length of the transversal confinement [3^. The latest 
is assumed to be very small such that the system is ef¬ 
fectively ID. The composite system where the only dif¬ 
ference between the atoms is the coupling strengths gij 
with respect to each other can be realized experimentally 
by using atoms of the same species in different hyperfine 
states. 

As shown in Fig. [T] we consider 3 atoms on a ring, 
with two of them being in the same hyperfine state and 
obeying fermionic or bosonic statistics, and the third 
(which we call the impurity) in a different hyperfine state. 
Therefore we have two important parameters, the cou¬ 
pling strength g between impurity and majority atoms 
and the g' for the bosonic case for the interaction be¬ 
tween the majority atoms. In the case where the major¬ 
ity atoms are fermions, the s-wave scattering is forbidden 
due to the statistics, and the only parameter that remains 
is the impurity-majority coupling strength g. 


Underlying Physics 

In Fig. [T] we schematically depict the extreme cases 
that indicate the basic underlying physical scenarios that 
appear in this system. Fig.[T](a) represents the integrable 
case where g = g' and the atoms are strongly repulsive, 
so they tend to maximize the distance between them. In 
Fig. [T] (b), only the majority atoms repel strongly each 
other (or obey the Pauli exclusion principle if they are 
fermions), and therefore they tend to localize on diamet¬ 
ric positions on the ring. The impurity is delocalised 
all over as a single particle since the interaction g with 
the others is vanishing. On the other hand, in Fig. [1] 
(c) the interaction between the bosonic majority atoms 
g' is vanishing or very weak, and the impurity strongly 
repels both. In this case the majority atoms tend to 
maximize the distance with the impurity, and this leads 
to an effectively “attractive” scenario, where they bunch 
together (this kind of effect has also been verified in har¬ 
monically trapped systems [1^). We will see that this 
indirect attraction of the majority atoms via the impu¬ 
rity induces various effects in the observable properties, 
the most striking being the peak in low momentum for 
the static structure factor. Let us stress here that at 
least these extreme situations are representative also for 
ensembles of many majority atoms, since the underlying 
physics are similar to those explained above. 

Ansatz 

Extreme cases in bosonic systems like the one depicted 
in Fig. [1] (a) allow also for exact solutions of the many- 
body Schrodinger equation that are obtained by means of 
Bose-Fermi mapping , which “maps” the infinitely re¬ 
pulsive bosons to free fermions. The Jastrow-type ansatz 
that we employ here is based on those solutions and also 
on the two-body integrable Lieb-Liniger case: 

— cos \hij (^\xi Xj\ 1/2)], (2) 

where kij (in units of 1/T) is related to the interaction 
strength gij by kij = 2 arctan(gy/2A:ij), with kij G [0,7r]. 
The Jastrow ansatz for our system of three bosons reads: 

Ip =JVC0S [k{\Xi -Xral \ “ 1/2)] 

X COS [k{\Xi - Xra2 \ “ 1/2)] 

X cos \k (jXmi Xm2l 1/2)]; (3) 

where A/" is a normalization constant, Xi denotes the po¬ 
sition of the impurity atom, Xmi and Xm 2 the positions 
of the majority atoms, and fc, k' are related to g and gP 
respectively. Although here we conveniently write the 
ansatz for our particular system it is clear that it can be 








3 



FIG. 2. (Color online) Comparison of results from Jastrow 
ansatz and exact solutions for different interaction strengths 
for (a) one-body correlation function and (b) pair correlations 
in the integrable bosonic case and (c) one-body correlation 
function for the impurity fermion and (d) pair correlation 
for a fermionic impurity-majority pair. Black dashed lines 
correspond to the results obtained using the corresponding 
exact wave functions. 


agreement with the existing exact solutions 3^. The 
same holds for the bosonic pair correlations as a function 
of the separation r = a;i — a; 2 , as shown in Figl2] (b). As 
expected, in the extreme cases of vanishing and infinite 
repulsion (Tonks-Girardeau gas) the agreement is at best 
since the Jastrow ansatz coincides with the exact solution 
[^ . For the case of two majority fermions with an im¬ 
purity, we compare the results of the Jastrow ansatz with 
the exact solution 32, 3^ in Figsl2] (c) and (d). Again 
we find that the agreement is very good in the weak and 
strong interaction limits, with slight deviations for in¬ 
termediate interaction strengths. Here p{x,x') refers to 
the impurity fermion, and the pair correlation function 
is calculated for an impurity-majority pair, so Xi = Xi 
and X 2 = Xml- For the homogeneous cases discussed 
here where an exact solution exists, the Jastrow ansatz 
provides a simple and reliable approximation of this so¬ 
lution. Therefore, we can state that even without a vari¬ 
ational approach, this ansatz captures the basic physics 
of the systems that we discuss here. We mention that an 
approach on correlations for bosons in a harmonic trap 
has been discussed in , while correlations for the two- 
component trapped fermionic system have been studied 


44, 45 


generalized for more atoms and different pairs. Indeed, 
for fermionic majority atoms, we modify the last term 
to sin [tt (xmi — Xm 2 )]^ which respects the fermionic ex¬ 
change property. Note that, contrary to the bosonic case, 
this term does not correspond to the exact ground state 
wave function for a pair of identical fermions (see ^ for 
details). However, as we will show next, this expression, 
besides of being simpler to handle, is able to reproduce 
the relative correlations for the cases that we study here. 

In Fig [2] we compare the results for correlations using 
the Jastrow ansatz to the exact ones given by the Bethe 
ansatz in the integrable cases. We particularly focus on 
basic quantities such as the one-body correlation func¬ 
tion. 


p{x,x')= J dX 2 ---dXN'4’ix,X2, ■■■,XN)'ip*{x',X 2 m--TXN), 

(4) 

and the pair correlation function 

g{xi,X2) = j dX3...dXN\lp{xi,X2,---,XN)\^, (5) 

where we assume that the A^-body wave functions are 
normalized to unity. In Figl2](a) we present the one-body 
correlation function with respect to the relative distance 
r = X — x' for the bosonic case for several interaction 
strengths. We hnd that the Jastrow ansatz is in good 


MEASURABLE QUANTITIES 
One-Body Correlation Function 


In our system we can find either the correlation func¬ 
tion for the impurity or the majority atoms by integra¬ 
ting out the coordinates of the other pair of atoms in each 
case. The one-body correlation function - and in partic¬ 
ular its off-diagonal terms - characterizes the degree of 
coherence of a given atom of the system. As we observe 
for the integrable case in Fig|3l(a) when the interactions 
are weak or intermediate (/c = 7r/2 corresponds to g = tt) 
there is coherence in the ensemble since diagonal and 
off-diagonal terms are almost equal. This behaviour re¬ 
sembles that of the Bose-Einstein condensate, which is a 
coherent ensemble of many bosonic atoms. On the other 
hand, when the repulsive interaction becomes very strong 
as in Fig [31(b) {k = Ott/IO corresponds to g « 35.7) the 
atoms tend to localise and lose coherence. The diagonal 
X = x' peak (and due to the periodic boundary conditions 
and symmetry the peak at distance L) becomes much 
more pronounced than the off-diagonal area around dis¬ 
tance L/2. This is typical for few body one-dimensional 
ensembles at strong interactions which are no longer rep¬ 
resented by a single coherent wave function but are oc¬ 
cupying higher single-particle states (orbitals), an effect 
known as depletion of the condensate or fragmentation 
3Si. We have a similar effect for the impurity atom 


when the impurity-majority interaction is strong and also 
for the majority atoms when the interaction between 
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FIG. 3. (Color online) One-body correlation function for representative cases interactions in the integrable, majority and 
impurity cases. Loss of coherence (depletion of the off-diagonal terms of p{x,x')) from (a) weak to (b) strong interaction in 
the integrable case, and from (c) weak to (d) strong repulsion from the impurity for the majority bosonic atoms. Behaviour of 
p(r = I® — x'l) as a function of the interaction parameter k in the cases of (e) integrable system, (f) impurity atom for strong 
majority interactions, majority atoms for (g) weak and (h) strong interactions and (i) majority fermions. 


them is strong. The most interesting non-integrable case 
is depicted in Figs [31(c), (d) where we see the loss of co¬ 
herence for the majority atoms solely due to the increase 
of the strength of the interaction with the impurity. Since 
the direct interaction term between the majority atoms 
is small {k = tt/IO, g w 0.1) one would expect that they 
remain coherent, but it is shown that if the impurity 
repels them strongly, then they get localised (see also 
schematic Fig[Tl(c)). Since the one-body density matrix 
is always symmetric with respect to the diagonal (as we 
see in Figs|3](a)-(d)) we can depict the whole behaviour 
as a function ot r = x — x' (as done in Fig|2]) and of the 
interaction. Therefore we show in Fig|3l(e)-(i) the whole 
behaviour of the correlation function as a function of k. 
The localization around r = 0 again shows that the loss 
of coherence appears gradually for high k. Very similar 
behaviour is obtained for the impurity atom as shown in 
FigO(f) due to the impurity-majority strong interaction. 
In Fig|3l(g) the bosonic majority atom is shown to lose 
coherence at large impurity-majority interactions even if 
majority-majority coupling is rather low. A comparison 
between the fermionic majority atoms and bosonic ma¬ 
jority atoms with infinite repulsion (k' = tt) shows that 
the fermionic atoms have more pronounced localisation. 
In other words, the infinitely-repulsive bosonic atoms are 
resembling fermions in their non-local correlation also as 
a function of the repulsive strength with the impurity, 
but they never get as radically localised as the majority 
fermions. 


Momentum Distribution and Contact 

An observable directly connected to the one-body cor¬ 
relation function is the momentum distribution: 

„(p) = /dxdxV(x,x>-«-'), (6) 

where p = 27rs and s is an integer. We show in the 
Appendix that this is equivalent to calculating first the 
many body wave function in momentum space for one 
atom 


(j){p,X2, ■■■,xn) = j dxie 'iI){xi,X2,...,Xn), 

and then integrating out the rest of the coordinates from 
its square modulus. 

It is known that the asymptotic behaviour (the large 
p tail) of the momentum distribution obeys an universal 
power law till, [ail: 


n(v) = —r, for p —>■ oo. 
P 


(7) 


An expression for C depending on the two-body correla¬ 
tion function is given in [H^: 


C = 


4(iV-l)p2(0,0) 


( 8 ) 


‘•ID 


where aiD is the one-dimensional scattering length, de¬ 
fined as ai£) = —2/mp, and p 2 ( 0 , 0 ) is the normalized 
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FIG. 4. (Color online) (a) Momentum distribution for differ¬ 
ent interactions, including integrable bosonic, majority and 
impurity bosons and fermionic cases (p is in units of 2n). 
Dashed black lines are the predicted values given by Eqs. [S] 
(integrable) and|9](non-integrable). (b) Contact as a function 
of k and k', obtained from Eq. [9]for the impurity atom, (c) 
Same as (b), but for a majority atom. 


two-body correlation function at vanishing distance be¬ 
tween the atoms (an analytical expression for p 2 in the 
case of three atoms is given in 37] , while thermodynamic 
limit formulas for weak and strong interactions are found 


51]). The proportionality constant C is called contact 


and is determined by the strength of the interaction. The 
contact characterizes all the universal properties of such 
systems that are independent of the details of the inter¬ 
action. For a fixed interaction, it is proportional to the 
probability that two atoms can be found at a short dis¬ 
tance from each other (in this case at zero separation). 
In Fig. HKa) we verify the Cjp^ asymptotic behaviour, 
not only for the integrable but for all cases considered 


here. We also verify that Eq. [8] is capturing the correct 
behaviour of the contact, for all our data points in the 
integrable case. For the non-integrable cases where inter¬ 
actions are different, there is no explicit formula in the 
literature for the contact. Therefore we write, in analogy 
to the calculations in (H^ . an analytical expression for 
the asymptotics of momentum distribution [52] : 



where we identify the pre-factor as the contact. The in¬ 
dex j is inserted to account for different interactions be¬ 
tween particle Xi and the other particles of the ensemble. 
For our particular case of three particles, the interaction 
between the impurity and the remaining pair is the same 
(so oJq = ajp), and the difference from the integrable 
contact will be carried only in the wave function. For 
the contact of the majority atoms, the scattering length 
will be different for each pair (af^ ^ Oid)- For all non- 
integrable cases this formula predicts very well the be¬ 
haviour of the momentum distribution at large large mo¬ 
mentum for our data points as is shown for exemplary 
cases in Fig. HJa). By using Eq. [5] we can then calcu¬ 
late the contact for all different interaction strengths, for 
impurity and majority atoms; the results for the contact 
as a function of the interaction parameters are presented 
in Fig. SDb) and (c). In general the value of the contact 
increases with the interaction strength, but the main ob¬ 
servation here is that the contact for the majority atoms 
is symmetrically dependent on k and k' while for the im¬ 
purity the major change comes, as expected, from the 
increase of the interaction with the majority atoms k. 


Static Structure Factor 

The structure factor is a property that defines how an 
ensemble of atoms scatters incident radiation. Exper¬ 
imentally, it is usually measured by two-photon Bragg 
scattering, in which the atom sample is subject to two 
detuned laser beams. The stimulated emission of light 
by the atoms gives rise to interference patterns that can 
be measured and contain information about the structure 
of the sample. For different wave vectors p and frequen¬ 
cies u! of the beams, it is possible to measure the dynamic 
structure factor S{p,uj) of the system. Its integral over 
all frequencies provides the static structure factor, de¬ 
fined for ID homogeneous systems as [li] : 

S'(p) = l-|- J dxj [g{xi,Xj) - 1], (10) 
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FIG. 5. (Color online) (a) Low momentum region of the static 
structure factor for several integrable and non-integrable 
bosonic cases, as well the fermionic majority-majority case 
(p is in units of 2tt). (b) Static structure factor as a function 
of k and k', obtained by using the majority-majority pair cor¬ 
relation function for bosons, with p fixed as 27r. (c) Same as 
in (b), but for an impurity-majority pair. 



where p is again quantised as 2tts and the pair correlation 
g{xi,Xj) is renormalized to iV(iV—l)/iV^ for convenience. 
For our system, the pair correlation function can refer to 
an impurity-majority pair (xi = Xi and Xj = Xmi) or a 
majority-majority pair (xi = Xmi and Xj = Xm 2 ), and 
N(N -1)/N^ = 6/9. 


In the low momentum region and infinitely repulsive 
regime, hydrodynamic theory 53| predicts that the static 
structure factor should behave linearly with p. Particu¬ 


larly, for the TG limit this behaviour is described by 
S(p) = |p|/27rn, where n = N/L is the particle den¬ 
sity. This result is characteristic of phonon excitations in 
many body systems, but we confirm it also in our few- 
body system in the integrable case as shown in Fig.jSKa). 
For high momentum the structure factor always con¬ 
verges to 1 (a high speed probe would not be scattered 
by the cold atom ensemble). For low momentum the be¬ 
haviour of the non-integrable cases deviates much from 
that of the integrable, the first peak being below the TG 
phononic behaviour or above it. 

Most importantly, for the majority bosons S can also 
be over 1, an effect that indicates a quasi-crystalline 
structure and is not present in the integrable case. It 
results specifically from the fact that the repulsion with 
the impurity is inducing an effective attraction between 
the majority atoms as depicted in Fig.[TDc). This is also 
a characteristic of the so-called Super-Tonks gas, which 
is an excited state on the attractive side of the Fesh- 
bach resonance that also exhibits a peak above 1 in the 
structure factor at low momentum as shown in (53| . We 
focus in Fig. [5Kb),(c) on the behaviour of the structure 
factor at low momentum p = 1. In (c) we observe again 
this peak arising smoothly for the majority atoms and 
being very pronounced at low k' and high k. For the 
impurity-majority pair, S' < 1, as is shown in Fig. [51(c), 
and the main change comes from impurity-majority re¬ 
pulsion. An important issue regarding the experimental 
determination of the structure factor is the possibility of 
measuring this quantity separately for impurity-majority 
and majority-majority pairs. This could be achieved, for 
instance, by having the beams polarized, in such a way 
that the photons could couple differently to the spin (or 
pseudospin) components in the system [55| . Another way 
to obtain this difference in the coupling involves a specific 
detuning of the laser beams, leading to separate measures 
in the components of the total structure factor, as argued 
in [H]. 


CONCLUSIONS AND OUTLOOK 

We have studied the non-local correlation functions 
of few-body bosonic and fermionic mixtures in a one¬ 
dimensional ring in the presence of an impurity in a dif¬ 
ferent hyperfine state. The interaction strengths between 
impurity-majority and majority-majority pairs may dif¬ 
fer, which renders the system non-integrable. The 
strength of this difference is responsible for pronounced 
effects in observable quantities. Most of the effects that 
we study on the one-body correlation function, on the 
contact and on the static structure factor refer to a 
weakly interacting majority pair which is strongly re¬ 
pelled by an impurity atom. In particular, this gives 
rise to effective attractive interactions and a pronounced 
peak in the structure factor. By means of our Jastrow 
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ansatz, which is in very good agreement with the ex¬ 
isting exact solutions for certain cases, we are able to 
provide results for several observable quantities in the 
non-integrable regime. 

The experimental verihcation of the quantities we 
study is possible since recent advances in ultracold atoms 
techniques allowed for measurements in few-body sys¬ 
tems. Many of our results also hold qualitatively for a 
larger number of majority atoms. Our approach can be 
extended to other cases which are relevant for experi¬ 
ments, like differently mixed systems and larger ensem¬ 
bles. 
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J 


The momentum distribution is usually obtained by performing the Fourier transform of the one-body correlation 
function p(x,x') (assuming a normalized wave function): 


i(p) = J dxdx'p{x,x')e 


(A.ll) 


We consider now p{x, x') as given in Eq. |4l By replacing it in Eq. lA.lli we get 


n{p) = J dx 2 dx 3 dx dx'tp{x,X 2 ,X 3 )'tp*{x',X 2 ,X 3 ) e 


(A.12) 


which can be rewritten as 


l{p) = J dX 2 dX 3 J dx'tlj{x, X 2 ,X 3 )e J dx' 'tp*{x' ,X 2 ,X 3 )l 


(A.13) 


From Eq.[71 J dxil:{x,X 2 ,X 3 )e = (j){p,X 2 ,X 3 ), so Eg. IA.13I results in 


n{p) = J dX2dX3\(j){p,X2,X3)\'^ 


(A.14) 


which can be calculated given that an expression for ^(p, 0 : 2 , 2 : 3 ) exists. Eor the momentum space representation of 
the wave function with respect to the coordinate xi in the bosonic case, this expression reads: 


y07^(j){p,X2,X3) 


' " "^2 + S3)]( - *(4e*P(i+2-=+-3)A:2 

_ 4 g®P(H- 3 : 2 -l- 2 a: 3)^2 _ g 2 ip{x 2 +X 3 ) _|_ ^ip{l+ 2 x 2 + 2 x 3 ) 

X COS [fc(a;2 — 2:3)] -I- /iie'‘P('^+^'^+^ 3 )_ ^ipx 3 '^j ^2 _(_ jjg)] 

_2gJp(l-l-a;2+a;3)(g*paD2 _|_ e’'P^^)kp{sill [k{x2 — 2:3)] d" sin [fc(l — X2 + 2:3)]) 
—ig'^v{x 2 +x 3 )p(^p(.Qg [/c (—1 -h X2 -I- X3)] — 2 jfcsin [fc(— 1 -I- X2 -b 2:3)])), 

__l_^g-zp(l+2.2+2x3) cos [fc'(i - 20:2 + 2x3)] ( - i(4e*P(i+2-3+"3)fc2 

_ 4 g*P(H- 3 : 2 -l- 2 a: 3)^2 _ g2ip(a:2+a:3) (_4^2 _(_ _|_ gip(l-|-2a:2+2a;3) 

X COS [fc(2;2 — X3)] — 4ig®p( 1+3^2+3:3)(g*pa:2 _ ^ipX 3'^^2 -\- X2 — 2:3)] 

_l_2g*p(H-3:2-i-x3)^gjpa::2 _|_ e'^P^^)kp{sin [k{x2 — 2:3)] — sin [fc(l -|- 3:2 — 2:3)]) 

, p{pcos [fc(— 1 + X2-\- 2:3)] — 2ifcsin [fc(—1 X2 + 2:3)])), 


if X2 > X 3 , 

(A.15) 


if X2 < X 3 , 


where Mb is the normalization for the bosonic wave function 


37| . For the fermionic wave function, we obtain 
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' 2(-4k^p+p3) e *^( 1 + 22 : 2 + 2 x 3 ) [7r(x2 - a;3)]( - k'^ 

_^gip(l+X2+2x3)^2 _ g2ip(x2+X3)^_^^2 -|_p2^ _|_ gip(l+2x2+2x3) ^_^^2 _|_p2^^ 

X COS [k{x 2 — 2:3)] + k{ 2 i{e^P^^ — e^P^^))kcos [fc(l — 2:2 + 2:3)] 

-{e^P^^ + e^P^^)p{sm [k{x 2 - 2:3)] + sin [fc(l - 2:2 + 2:3)]) - ie^ip(^^+=^3) 


3/N^(I){P,X2,X3) = < 


x(-l + e"P)p{pcos [fc(-l + X 2 + 2:3)] - 2iA:sin [A:(-l + 2;2+x3)])), 

_^__i_^g-zp(i+2x2+2x3) sin[7r(x2 - 2 : 3 )] ( - i(-4e*P(^+2“2+x3);i,2 
_|_4gip(l+a:2+2x3) ^2 _ g2ip(x2+X3) _|_p2^ _|_ gip(l+2x2+2x3) j'_^^2 -|-p2^^ 

X COS [k{x2 — 2:3)] — cos [k{l + X2+ x^)] + 

x(g*pa :2 _j_ e*P“^))fcp(sin [A:(a;2 — 2:3)] — sin [fc(l + X 2 — 2^3)]) — 

. x(-l + e^P)p{pcos [fc(-l + X 2 + X 3 )] - 2ifcsin [k{-l + 2 ; 2 +x 3 )])) 


if X2 > X 3 , 

(A.16) 


if X 2 < X 3 . 


The normalization of the wave function for the fermionic system is given by 

2fc® — 2fc^7r^ — 4fc^7r^ + + 2k^Tr'^ + 2k^TT^ cos (2k) — tt^ cos (2k) + 4 /c® sin (fc) — Sk^ir^ sin (fc) + 4fc7r^ sin (fc) 


A/f = 


16A:^(—fc + 7r)^(fc + tt)^ 







